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In this paper the effects of vibrations at high frequencies onto a freely falling two-body system in 
Schwarzschild spacetime are investigated. As reference motion of the same system without vibrations 
a circular orbit around the central body is considered. The vibrations induce a perturbation on this 
motion, whose period is close to the orbital period, in agreement with the simpler situation of the 
Shirokov effect [J]. In general relativity the amplitude of the perturbation is dominated by high 
■ velocity effects, which grow linearly in the radius r of the circular orbit, while the leading term 

' surviving the Newtonian limit decays as 1/r. Thus even for very large radii a significant difference 

between Newtonian physics and general relativity is found. We give an estimate of this effect for 
• some molecular vibrations of a system orbiting around the Earth. 
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O ! I. INTRODUCTION 

CT 

Changing the trajectory of an orbiting spacecraft can be done with internal motions. This is due to the work of the 
' tidal forces experienced by the system. Two-body systems orbiting the Earth have been studied in [1, 0, 0] . As most 
important limitation, large effects in Newtonian physics are resonant effects, where the frequency of change of shape 
is linked to the motion of the system (e.g. its orbital frequency or its orientation.). More recently Wisdom Q and 
Gueron et al. [y, LZ[ studied similar situations within general relativity. Within this more general setup non-resonant 
effects can become large. Thus, relevant deviations from the free motion can be obtained for frequencies of the change 
\ of shape completely independent of the orbital motion. 

In this paper we investigate a vibrating two body system similar to the one of Ref. [7[. As reference motion 
a circular orbit is considered (as opposed to a radial free fall as in Ref. 0) and it is shown that this leads to a 
systematic deviation between the Newtonian result and general relativity even for large radii of the orbit (i.e. weak 
gravitation) since the leading perturbation in general relativity is a high- velocity effect, where the vibrational velocity 
acts as the key parameter. In particular, the maximal deviation per orbit is found to increase with increasing radius 
of the orbit. Rather than with the result of Ref. [?J our result is related to the so-called Shirokov effect [l|. Still, the 
. ! latter work just considers one perturbation at a single time instance and not a continuous perturbation as done in 
this work. Accordingly we find effects which are at least three orders of magnitude bigger than the effect reported in 

a. 

The paper is organized as follows: the following section presents the model and the geometry used throughout. In 
Section IIIII a simplified model is studied in detail and it is shown that in general relativity the total deviation from 
the free motion increases linearly with the distance from the central body whereas the dominant term surviving in the 
Newtonian limit decreases. In the subsequent section, we make a numerical study of the full system and we show that 
the effect found for the simplified case still dominates the behavior. In Section [V] we comment possible experiments 
and finally our conclusions are presented in Section [VTl 
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FIG. 2: Profile of the constraint. 



II. THE MODEL 

Following the previous works @, 0, 0|; & vibrating or oscillating system is implemented as a collection of point 
masses whose relative positions are related by time dependent constraints. The specific model used here essentially 
is equivalent to the one of Ref. 7]: a two-body system made from two test masses connected by a massless tether, 
whose length l(t) is imposed by an oscillating constraint. To further simplify the situation the two masses are always 
considered to be equal and put to unity in all calculations. This system shall orbit around a central body, which is 
described by means of Schwarzschild geometry. We assume that the two masses have the same orbital plane, where 
the oscillations take place; therefore in Schwarzschild coordinates, the motion of each of the two point masses can 
be defined in terms of the radial and azimuthal coordinates (r^ipi) (where i — 1,2), while the polar angle angle is 
dropped. The system may be fully described in terms of the variables {r\, tpi), the relative angle 9 — ip2 — <pi and 
the constraint l(t) (see Fig. (Ha)). However, in many situations it is useful to use instead the coordinates of the 
geometrical center of mass (r, ip), the relative angle between r and I, (3, and the constraint l(t). This situation is 
depicted in Fig.QJb). 

In Ref. Q it has been shown for a radial free fall that asymmetrical oscillations induce a quite different change in 
the motion of this system than symmetrical do. To assess the same question for closed orbits, a similar profile of the 
constraint including a potential asymmetry is used here; however, it will be shown below that the exact profile of the 
constraint is unimportant. Our constraint is described by four parameters: its frequency Cj = its amplitude 5i, its 
minimum length Iq and an asymmetry parameter a. The asymmetry parameter a, taking values in the range [—1, 1], 
indicates how much the constraint fails to be symmetric with a = being the symmetric casepj]]. Unlike Ref. 0, we 
need to introduce a minimal length Iq between the two masses in order to avoid divergent angular momentum of the 
spinning two-body system. The profile of the constraint is shown in Figure [2] 

The constraint I relates the coordinates of the two masses as 

I 2 = r\ +r\ - 2rir 2 cos6» . (1) 
Since l/r and consequently also 8 is always kept small, r 2 in very good approximation can be expressed as 

ra = n + ^l 2 (t)-r\0^ = n + Sr 12 . (2) 

In many explicit calculations the above system is yet too complicated. To reduce clutter a simplified situation will 
be considered. In this system an additional constraint is imposed which enforces ipi = tf2 = l P and thus 8 = 0. We 
call this system the double constrained system (by opposition to the full system). 

We point out that no corrections from spacetime curvature have been taken into account in the calculations of the 
relative coordinates between the two masses. Though our ansatz mainly will be justified by the final result, we would 
like to comment about these issues at this place already. Most importantly it should be realized that we intend to 
place the constrained system exclusively in regions of weak gravity, where r s /r <C 1, and thus the ensuing errors are 
expected to be small. Furthermore it will be shown in Section IIII Bl by means of the expansion of the equation of 
motion of our vibrating system that a systematic difference exist between general relativity and the Newtonian theory, 
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FIG. 1: Representation of our model and the different vari- 
ables used. 
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which is not affected at all by this simplification. Thus, this small error just affects the numerical results presented 
in the paper; our main conclusion based on analytical expansions remains unchanged. 

Still, the interpretation of the implementation of the constraint and the resulting deviation Ar from the reference 
motion in terms of the Schwarzschild coordinate r may be questioned. However, we are not interested here on 
finding an experimental protocol to implement a particular constraint, or to describe a realistic constraint. If so, 
one could implement the constraint in the Fermi- Walker coordinates [1[ of the reference motion, which in good 
approximation are comoving coordinates of the vibrating system. This would not change our results qualitatively. 
Moreover, the deviation Ar is a coordinate distance between two defined events along the perturbed and the non- 
perturbed trajectories. We will express it in term of radar distance, which is then independent of the particular 
coordinates used in the calculations. 

Let us now write down explicitly the relevant actions of our system in general relativity. As explained above, we 
assume the spacetime geometry to be described by the Schwarzschild metric, and we implement the constraint in 
Schwarzschild coordinates. Thus the action becomes fl2l: 

S = - [ dt 



ii + V^2 + A(r2-ri-(5ri 2 ) , (3) 



where A is a Lagrange multiplier, t is the Schwarzschild time and 



I -) (l -) r: rf^;. (4) 



r, J \ r. 



with r s being the Schwarzschild radius of the central mass, and () = d/dt. 

In the following sections, we study the effect of the oscillations of the length of the constraint, l(t), onto the motion 
of the two-body system. We compare the trajectory of the vibrating system with a reference motion, which is the 
motion of the same two-body system with identical initial conditions but without vibration. For the reference motion 
the simplest situation, namely circular orbits, are chosen and the two-body system always is considered to be aligned 
in the radial direction. We emphasize that the reference motion is not a geodesic. However, it is straightforward to 
show that 4> is constant for the reference motion. The initial conditions of the non- vibrating and the vibrating systems 
are denoted by: r(0) = r , <p(Q) = 0, 9(0) = 0, f(0) = 0, ip(0) = uo and 9(0) — 0. By imposing the reference motion 
to be circular we get from the conservation of angular momentum a relation between ro and w; for the following the 
explicit relation between tq and u> is not important so we do not reproduce it here. 



III. THE DOUBLE-CONSTRAINED SYSTEM 



We first study the double constrained system, which is described by the action ([3]) with the additional constraint 
9 = 0. It turns out that this system is much simpler than the full system and even allows to make analytical 
expansions, since 

• it contains one degree of freedom less than the full system; 

• it allows a much simpler implementation of the constraint l(t); 

• it includes only two frequencies (the frequency of the constraint l(t) and the orbital frequency) instead of three 
in the full system, which in addition includes the frequency of the oscillation of the two-body system around its 
geometrical center of mass. 

We hasten to add that the motivation to study the double-constrained system is of purely theoretical nature, mainly 
the fact that this system can be treated analytically. We do not claim at this point that this specific system easily 
can be implemented in a real experiment. Nonetheless, it will be shown in Section ITVl that the differences between the 
double-constrained and the full system remain very small. This justifies to study the double-constrained system more 
in detail since the dynamics of the full system are clearly dominated by the same effects as the ones to be worked out 
in this section. 



A. Linearization of the equations of motion 



Of course, the exact equations of motion can be integrated numerically, which however can be quite time consuming, 
especially for large radii. Without use of extensive computer resources it is not possible to integrate the equations 
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for a system orbiting the Earth. Fortunately, the double constrained system can be treated analytically in very good 
approximation. To do that, the equations of motion are linearized around the reference motion, which we have defined 
to be the motion of a non-vibrating system with the same initial conditions. As pointed out above, this reference 
motion is assumed to be a circular orbit characterized by the radius ro and the orbital frequency uj. The parameters 
of the vibrating system become r = ro + Sr and tp = cut + Sip and yield equations in the deviation variables Sr and 
Sip, which are then expanded to first order. One of the two resulting equations is of the type 



A{t)Sr + B{t)5(p = d (t)Sr + D x (t)Sr + E l (t)5<p + F 1 (t) . (5) 
Since the angular momentum 



is conserved, the linearized angular momentum allows to substitute Sip in (O, which yields 

A(t)Sf + B(t)5Cp = C{t)Sr - D(t)5r - F(t) . (7) 

Here A(t), B(t), C(t), D(t) and F(t) are complicated functions but they are all periodic with the period of the 
constraint T = 1/u). At this point it is important that our system only includes two frequencies, w and the orbital 
frequency uj. As is confirmed by integrating the exact equations, Sr oscillates with a period of the order of uj. Since 
ui « ui, Sr does not relevantly change during one period T and thus all coefficients A(t), B(t), C(t), D{t) and 

F{t) may be averaged over one period of the constraint, i.e. A(t) is replaced by A — ~ A(t)dt etc. Using this 
simplification, one can show that B = C — 0, and the solution of the equation of motion simply becomes 

p 

Sr = ^ (cosuj't - 1) , (8) 

where uj' 2 = D/A. A numerical integration shows that D/A is positive, and that ui' ~ u. Moreover, \Sr\ <C ro implies 
that F/(Dro) <C 1. Then, the resulting trajectory of the vibrating system is an ellipse of eccentricity e w F/(Dr ). 
The advantage of this formulation is that we only need to perform an integration of the three functions A, D and F 
over one period of oscillation instead of integrating the full equation ([7]) over the desired time of evolution, e.g. one 
revolution around the central body. In this way we directly obtain an expression for the maximal value Ar of the 
deviation from the reference motion: 

IF 

Ar = ^ . (9) 
D ' 

The corresponding radar distance Al can be obtained by integrating the relation dl 2 = (1 — rs/r) _1 dr 2 between ro 
and ro + Ar. The radar distance does not depend on the particular coordinate system we used for the calculation, 
and in principle can be determined by instruments on-board the spacecraft. 



B. Analytical calculation 

The result still does not allow a full analytic treatment, since the functions F and D cannot be integrated 
analytically. Still, further expansions of the equation of motion ([7]) are possible since the two parameters l(t)/ro and 
l/c both remain small [l3j . We do not expand in the parameter r s /r as this is not necessary to obtain our result. 
Notice however that all our discussions take place in a regime of weak gravity, where r s /ro -C 1. 

In our calculations, l/c typically is of the order of 1CP 4 -10~ 5 , while //ro ranges from 10 -5 at relatively small radii 
up to 10~ 12 in Earth orbit (cf. also the numbers presented in Sect.fV}. Then we obtain: 

Ar _ S 2 Cj 2 Si(a) 



c2 (^X 1 -^) 



5M ^ f 2^S 2 (a) + *S 3 (a) ) (10, 



ro r 



+ O 




FIG. 3: Numerical evaluation of Si (a), 82(a) and S3 (a). 



Si, S2 and S3 are functions that depend on the asymmetry parameter and are represented in Fig. [31 The deviation 
Ar contains two kinds of terms: a high velocity effect on the first line and a purely gravitational effect on the 
second line. These two effects behave differently with ro- In the following, we will assume that l/r s < 10~ 3 and 
r s /r <C 0.01 — 0.1 [l4j]; then the high velocity effect is much bigger than the purely gravitational effect. This 
dominance increases with increasing radius, since the velocity effect grows linearly with tq (la |. while the purely 
gravitational effect decays like l/r . Keeping only the leading term in ITU)) while expanding in r s /r , we obtain: 



The corresponding radar distance, to zeroth order in r s /ro, is A^ ~ Ar. Then, to this order the deviation (fTTj) does 
not depend on the choice of the coordinate system we used for the calculation. Also it is seen that the deviation is 
determined solely in terms of the change of the length of the constraint, Si, and its frequency; as mentioned already 
above the exact form of the constraint remains unimportant. 

The Newtonian limit of equation (flO)) is obtained by taking the limit l/c — > and r s /ro — » 0: 



Given the fact that the leading effect in general relativity is a high velocity effect it is evident that the differences 
between general relativity and Newtonian gravity are much more fundamental than their mere difference in the 
interpretation of the radial coordinate (we point out again that all the results apply in regions with weak gravitation 
and small velocities of the system as a whole, where this difference is small). Since the effect due to the velocity 
term grows with increasing radius, while the difference in the interpretation of the radial coordinate decreases, this 
conclusion applies to any desired accuracy if only gravity is weak enough. 

Before studying Ar in detail, we want to show that the resulting motion of the vibrating system to a very good 
approximation is an ellipse. From ([8]) the frequency or period of the perturbation Sr can be calculated. With T being 
the time of revolution of the reference motion, the period of the perturbation is found as 



where ki is a function of the ratio r s /ro. Note that ki — » 1 when r s /ro — > 0, so the perturbation in the orbital period 
stays even in a weak gravity field. To the same order in the expansion, the advance of the periapsis is given by 




(12) 




(13) 




(14) 
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FIG. 4: Ar as a function of ro for lo = 0.05, a = 0, Si = 5 x 10 -3 and lo = 5x 10 -3 : comparison of the numerical simulations 
and the analytical expansion (|10l ) In addition the leading purely gravitation term has been included in the figure to illustrate 
the difference between general relativity and Newtonian gravity. 

where hi is again a function of r s /ro, which for large radii is of order 1. At large radius, the second term dominates 
over the first one but in absolute terms stays very small (cf. Sect.|V]for further discussions). Thus the resulting motion 
indeed is an ellipse with a small advance of the periapsis. 

It is interesting to look at the conservation of energy. In general relativity the Hamiltonian H = Hi + H2, with 

Hi = ^-==t, (15) 

V -^i 

may be used as a definition of "energy" . Obviously, this Hamiltonian is conserved for the non- vibrating system. 
Together with oscillations, the energy changes during each cycle of the oscillation. Nonetheless, looking only at the 
starting or end points of the vibrations, t = (0, T, 2T, . . .), the energy remains conserved to second order in l/rg and 
to second order in l/c. This means that the variation of energy (if there is any at all) is of higher order and therefore 
remains very small. This analysis is confirmed by numerical simulations; however, due to numerical errors, it is hard 
to quantify its actual change. 

In conclusion the vibrations do not change the energy of the system considerably. Thus, although the orbit of the 
vibrating system geometrically is an ellipse, the motion of the system is different from the geodesic motion described 
by the same ellipse, since the energies (and thus the velocities) of the two situations are different. 

In the remainder of this section we comment on the relation of our system to the so-called Shirokov effect 0, [t| [l(| • 
In the Shirokov experiment, a non- vibrating and point-like body is perturbed from a reference circular orbit. Then 
its trajectory oscillates around its reference orbit, with a different period if the movement is along the radial or the 
zenithal directions. The experiment consist to find the difference between these two periods. For a radial perturbation, 
Vladimorov has shown that the oscillation is due to the quasielliptical nature of the orbit, with a general relativistic 
Mercury type displacement of the periapsis [lCj . The period and the advance of the periapsis of this orbit can be 
compared to the one we find for the vibrating system, except that for the latter there are supplementary corrections 
due to the oscillation velocity. However, these two orbits are located differently compared to their circular reference 
orbit: the orbit of the non-vibrating body in the Shirokov experiment intersects with reference orbit twice per orbital 
period, whereas the orbit of the vibrating system crosses its reference orbit only once per orbital period. This is due 
to the small correction —F/D contained in Equation (jSj). 

Moreover, in [lj, Shirokov considers as given the maximum deviation of a trajectory perturbed in the radial direction 
of a non- vibrating body from its reference circular orbit; this is very different from the situation we describe where the 
maximum deviation is the quantity we want to determine. In fact, in the Shirokov experiment, the only thing that 
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FIG. 5: Comparison of the numerical simulations with the analytical result : Ar as a function of ui for different values of 
a — {0,0.15,0.4}, ro = 200, Si — 5 x 10 -3 and lo = 5 x 10~ 3 . The results for negative a are exactly the same as those for 
positive values. 



differs between the reference and the perturbed body are their initial conditions, whereas in our situation the initial 
conditions of the non-vibrating and the vibrating bodies are the same; the difference is in their internal movement. 
Then Shirokov calculates the coordinate distance £ between the body perturbed in the radial direction and a body 
perturbed in the zenithal direction after n orbital periods. For n = 10 orbits he obtains £ ~ 10~ 6 cm. This is at least 
3 orders of magnitude below the effect we find for vibrating molecules after one orbit (see section V). 



C. Comparison and results 



In this section results from the two presented approaches, the numerical integration of the complete equations of 
motion and the analytical result from a systematic expansion (|10p . shall be compared. 

Figured] presents [16[ in a log-log scale the numerical analysis of the full equations and the analytical result for Ar 
as a function of ro . For comparison the Newtonian result is included as well. As can be seen, the general relativistic 
result is perfectly linear above ro ~ 100; below that value strong gravity effects appear due to the corrections in r s /ro 
in the relation IjlOp . Furthermore it is seen that the analytical result fits perfectly well the numerical data. Considering 
the Newtonian result the curve is linear in a log- log scale which fits the 1/ro behavior as found in Eq. (fT!Z]) . We can 
clearly see on this figure the divergence between the classical approach and the relativistic one. 

In Figure [5] the variations of the frequency of the constraint, uj 1 and the asymmetry parameter, a, for fixed radius 
ro are presented. Again the analytical result is compared with the full numerical integration and the quadratic 
behavior in the frequency is confirmed. As can be inferred from the evaluation of S% (a) in Fig. [3J Ar increases if 
the asymmetry parameter is chosen different from zero. However, the figure clearly displays that an increase in the 
frequency is preferable to increase the effect. It is important to point out the differences of this result compared to 
a radial fall of the same system Q: there, the asymmetry parameter plays a central role since no effect is found at 
a = 0, furthermore the result in that situation is asymmetric for a — » —a while it is symmetric here. In conclusion 
one can say that the potential for asymmetric constraints in circular orbits is very limited. We also mention that 
the authors of Ref. Q found a characteristic plateau when varying the frequency. It should be pointed out that this 
"plateau" in the current context should rather be seen as a linear behavior since here we integrate the effect over a 
certain fixed time rather than a fixed number of oscillations as done in Q • Still, the circular motion shows a stronger 
dependence on the frequency, since it increases like uj 2 for increasing frequency. Indeed, the effect found in Ref. [7J 
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FIG. 6: Comparison of the numerical simulations with the analytical expansion (|10[) : (left) Ar as a function of Si for ro = 200, 
Q = 0.05, a = and l = 5 x 10 -3 ; (right) Ar as a function of l for r = 200, Cu = 0.05, a = and Si = 5 x 10~ 3 . 



is of order O(lro/c 2 ) oc tD, where ro is the velocity of the radial free fall of the non-vibrating system, to be compared 
with the effect from the circular case, of order 0(l 2 /c 2 ) oc Co 2 . 

Also, we can compare the results from varying the amplitude of oscillations Si. Figure UlI CI shows Ar as a function 
of 5i obtained from numerical simulations and from the analytical expansion. Again the analytical result fits perfectly 
well the numerical simulation, confirming the quadratic behavior in this variable 

Finally, the behavior of Ar for different values of is presented. Figure UlI CI reproduces the linear behavior in this 
variable with a good agreement of numerical data and analytical expansion. Comparing the Figures Till CI and IIII CI 
it is seen that has a small influence on the result in comparison to 5i. Thus in practice it is important that the 
system vibrates with a large amplitude, while the minimal distance of the two point masses is of minor relevance. 



IV. THE FULL SYSTEM 



After a rather detailed analysis of the double constrained system, some results about the full system are presented. 
In this section the constraint 9 = is relaxed. In other words the two-body system needs no longer be aligned radially, 
but can freely spin. Still, as initial conditions a radial alignment is chosen. With this choice of initial conditions the 
angle (3 (see Fig. [T]) will oscillate at a characteristic frequency. It is no longer possible to get analytical formulas as in 
the previous section, since in the equations of motion all variables are coupled and it is no longer possible to decouple 
them. Still, we can integrate numerically for rather small radii the full equations of motion. 

In these simulations the effect seen in the double constrained system, namely the difference in the radius Ar, is 
found as well. This result is presented in Fig. IIVI where the solid line is the analytical expression for the double 




FIG. 7: Comparison between the numerical simulations for the full system and the analytical expansion (|10[) for the double 
constrained system: (left) Ar as a function of ro for Co — 0.05, a = 0, Si — 5 x 10 -3 and io = 5x 10 -3 ; (right) Ar as a function 
of Co for different ro with a = 0, Si = 5 x 10~ 3 and lo = 5x 10~ 3 . 
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r (GM/c 2 ) ffl (c 3 /GM) 

FIG. 8: Representation of A/3 and uip for the full system : (left) as a function of ro for uj — 0.05, a — 0, Si = 5 x 10~ 3 and 
Zo = 5 X 10~ 3 ; (right) as a function of u> for different ro with u> = 0.05, a = 0, Si — 5 x 10~ 3 and io = 5x 10~ 3 . 

constrained system. As can be seen, the difference between the two systems remains unimportant. The same applies 
if Ar is shown as a function of the frequency of the constraint Co. 

This result is shown in Figure IIVI and again the analytical result from the double constrained system fits perfectly 
well the simulation points of the full system. Thus, as far as Ar is concerned, the above results of the double 
constrained system describe this more complicated situation very well. 

Still, on top of the shift in the radius the relative position of the system, described by the angle /3, will oscillate 
around its equilibrium position /3 = 0. In the following, we are mainly interested in this new oscillation, which is 
characterized by its amplitude A/3 and its frequency cop. 

Figure HVl shows the behavior of A/3 and cop with the radius r . They are decreasing with the distance from the 
central body, as can be expected on general grounds since the gravity gradient decreases with increasing radius. As 
a consequence, we expect that at large radii the full system will behave like the double constrained system. 

Finally, Figure [TVl shows the behavior of Ar, A/3 and top as a function of the oscillation frequency Q. It is seen 
that the amplitude of oscillation of the system around its center decreases with Co but the frequency of oscillation cop 
seems to increase linearly. 

V. SOME COMMENTS ON POSSIBLE EXPERIMENTS 

Could the effect described in this paper be measured on-board a spacecraft orbiting the Earth? Since we mainly 
deal with a high velocity effect, the weak gravitational field does not form an obstacle. On the contrary, the maximum 
deviation from the reference motion grows linearly with the radius, which on the other hand is at the price of an 
increased integration time (half a revolution period). A promising system to measure the effect could be molecule 
vibrations; in the following we estimate the result from such a system for an idealized situation, where a single 
molecule is placed in orbit with its vibration direction exactly aligned radially. Using a semi-classical model, we can 
describe a molecule as a quantified anharmonic oscillator. With this simple model, we find the values of the oscillation 
parameters as presented in Table [I] 
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TABLE I: Value of the oscillation parameters for two different molecule vibrations. M is the mass of the Earth. 

Notice that these values are given in geometrical units, with M the mass of the Earth; therefore specific numbers 
are only valid in this environment. Nevertheless, the leading contribution in (jTUJ) only depends on the velocity SiCo, 
which is independent of M. Therefore the results presented here to a very high accuracy hold in a (weak) gravitational 
field different from the one of the Earth. 
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FIG. 9: Ar as a function of the initial radius rn for the double constrained system with values corresponding to the ground 
state molecular vibration of HBr and CO (see Tab|T]) For the simulation of the linearized system Q a central body whose 
mass is equal to the Earth's mass (M = M©) has been assumed, the analytical result is independent of M. 



If the molecule is taken in its ground state, the result as depicted in Figure [5] is obtained, where both, the analytical 
expression as well as simulation points for the linearized system, are plotted. We can see that for Earth radius and 
low Earth orbits (r ~ 1.43 10 9 GM/c 2 ), the effect is of the order of 

Ar - lQ- 2 GM/c 2 - 5 x 10~ 2 mm. (16) 

The result does not look promising, but there exist different strategies to increase the effect. Firstly, one might try to 
work with excited states, which can increase the effect by about a factor of 100. Secondly, the experiment could be 
placed in a higher orbit. As an example, one can win a factor of 10 by placing the experiment in a geostationary orbit. 
Taking both strategies together, an effect of about 5 cm results. Of course, one might also try to place a spacecraft 
in an orbit around a central body different from the Earth. If the experiment is placed in an orbit around the sun 
with a radius equal to the one of the semimajor axis of the Earth's orbit, we get a displacement of about 10 m for 
the ground-state vibrations, and about 1 km for excited states. 

As has been mentioned in Sect. IIII Bl besides the shift Ar the perturbed trajectory is characterized by an advance 
of the periapsis of the ellipse. We mention that for molecular vibrations the typical values for this advance are of the 
order of 10~ 12 rad and thus not of experimental interest. 



VI. CONCLUSIONS 



In this paper we investigated the general relativistic effects of vibrations on a two-body system placed in an orbit 
around a central body. In our analysis the background spacetime was taken as Schwarzschild spacetime and the 
reference motion of the two-body system was assumed to be a circular orbit. Within a suitable expansion of the true 
motion around the reference motion, the equations of motion can be solved analytically, which led to our main result 
the relation (fTU|) . It was found that the vibrations deform the circular orbit to an ellipse, described by the maximal 
deviation of the reference motion and an advance of the periapsis. We stressed the similarities and differences between 
the experiment we described with the one by Shirokov [l[. 

In general relativity, the dominant contribution to the deformation is a high-velocity effect and the maximal devi- 
ation from the reference motion grows linearly with the radius rn thereof. In contrast to this, no high-velocity effects 
are present in the Newtonian theory and the dominant effect decays like 1/rn. Since the relativistic effects are at large 
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radii orders of magnitude larger than the Newtonian one, it should be possible to measure this effect in the Solar 
system. 

An estimate of the orders of magnitude has been presented for molecular vibrations. A promising strategy could 
be to place a spacecraft in a high orbit (geostationary or higher) and to use excited states as vibrations instead of the 
ground state. Of course, many questions considering possible experiments are open. Here we considered an idealized 
situation where a single molecule is placed in a circular orbit. Moreover, we did not take into account that the two 
atoms within the molecule can have a different mass. More realistic situations, possibly including a statistical analysis 
over a large number of vibrating molecules, should be analyzed. Also, different experimental setups could be studied, 
e.g. vibrations in crystals or nanoparticles. However, for most of these situations more complex models than presented 
here should be studied. Further topics could be the use of more complicated background spacetimes and different 
orbits as reference motions. 
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